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ABSTRACT 

The  theoretical  solution  for  the  flow  beneath  V/STOL  aircraft  was 
extended  to  include  tilted  jet  configurations.   A  laboratory  model  was 
constructed  to  test  the  effect  of  variation  of  the  parameters  governing 
the  flow.   Free  streamline  plots,  pressure  coefficients  on  the  ground 
and  fuselage  and  velocity  profiles  in  the  nozzles  were  determined  from 
hot-wire  anemometer  traverses  and  micromanometer  readings.   Experimental 
data  compared  favorably  with  the  theoretical  determinations. 
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SYMBOLS 

A  a  constant 

C  pressure  coefficient 

P  F 

I  flow  region  in  the  Q    plane 

K  complete  elliptic  integral  of  the  first  kind 

K'  defined  by  K' (k)  =   K(k') 

k  modulus    of  elliptic    integral 


h 


k'  defined   as   \'l   -   k 

p  pressure 

p  ambient    pressure   outside   jets 

Q  an  angle   used   to  evaluate   Q  functions 

R  Reynolds    number,    AV   /v 
e  o 

T  complex  variable    in   T-plane,    §  +   ir] 

t  complex  variable    in    t-plane 

u  velocity   in  x-direction 

V  velocity 

VT  velocity  at   nozzle   exit 

V  average   velo   ity  at   nozzle   exit 

V  velocity  along   free    streamlines 
o  J  & 

v  velocity   in   y-direction 

W  complex   potential,    §  +   iijr 

x,y  coordinates    in   physical   plane 

z  complex  variable,   x  +   iy 

a  tilt   angle    of   nozzle 

3  plate   angle 

T  flow  region  in  the  T-plane 


A  width  of  nozzle 

6T,  6D  asymptotic  widths  of  streams  flowing  to  left  and  to  right 
L   R 

Q  complex  conjugate  velocity,  u-iv 

9  Theta-function 

®  function  defined  in  equation  (4) 

§,r|  coordinates  of  T-plane 

p  density 

$  velocity   potential 

i|i  stream   function 

fi  function  defined    in  Appendix  A 

Subscripts : 

A,G,H  refer  to  stagnation  points 

r  refers  to  boundary  point 
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I .   INTRODl  CION 

In  recent  years  crowded  airport  facilities,  the  search  for  rapid 
urban  transportation  and  the  need  for  military  aircraft  capable  of 
operating  from  advanced  bases  has  resulted  in  a  great  deal  of  attention 
being  given  to  V/STOL  (vertical/short  take-off  and  landing)  aircraft. 
Much  research  has  been  performed  on  various  models  and  full-scale  versions 
of  these  aircraft.   Some  operational  types  have  been  developed. 

When  operating  near  the  ground  these  aircraft  are  quite  similar  in 
principle  to  ground  effect  machines.   Ground  effect  can  be  either  beneficial 
or  detrimental  depending  upon  the  configuration  of  the  aircraft.   Land 
erosion,  damage  to  objects  near  the  operating  area  and  recirculation  of 
debris  and  hot  exhaust  gases  beneath  the  fuselage  are  some  of  the  prob- 
lems involved  with  hovering  near  the  ground.   The  reingestion  of  exhaust 
gases  reduces  the  efficiency  of  the  engines.   Recirculating  debris  may 
also  enter  the  engines  or  severely  damage  the  underside  of  the  aircraft. 
Some  pilots  have  experienced  handling  difficulties,  drastic  loss  of  power 
and  visibility  problems  as  a  result  of  ground  effect. 

The  recirculating  flow  underneath  the  aircraft,  however,  can  also 
produce  a  lifting  force  which  can  significantly  reduce  the  thrust  necessary 
to  operate  the  aircraft.   Small  modifications  to  the  aircraft  such  as 
tilting  of  engines  can  often  reduce  the  unfavorable  effects  to  such  an 
extent  that  the  additional  lift  provided  by  the  recirculating  flow  can 
produce  improved  performance.   A  mathematical  model  of  the  flow  field 
beneath  V/STOL  aircraft  would  permit  designers  to  estimate  the  effects 
of  various  modifications  analytically. 


Unfavorable  ground  effect 


Favorable  ground  effect 


Surface  erosion  and  debris  recirculation 
FIGURE  1 
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II.   OBJECT 

The  theoretical  solution  for  the  flow  fie  Id  beneath  V/STOL  aircraft 
with  two  vertical  two-dimensional  jets  has  been  presented' in  a  series  of 
three  NASA  reports.   The  purpose  of  this  research  was  to  obtain  experi- 
mental verification  of  the  present  theory  and  to  extend  the  theory  to 
include  other  possible  configurations. 
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III.   THEORY 

An  analytical  solution  to  the  flow  beneath  a  V/STOL  aircraft  was 
obtained  by  Goldstein  and  Siegel  [Ref  3].   They  considered  a  two- 
dimensional  model  of  an  aircraft  with  a  slot  jet  exhaust  or  fan  pod  in 
each  wing.   The  analysis  was  performed  assuming  an  incompressible, 
inviscid,  isothermal,  irrotational ,  steady  flow.   A  flat  plate  at  an 
angle  to  the  ground  represented  the  fuselage.   The  flow  field  was  mapped 
into  a  hodograph  plane  and  a  complex  potential  plane.   A  brief  outline  of 
this  analysis  is  given  in  Appendix  A. 

So  that  entrainment  could  be  ignored  this  solution  assumed  that  the 
nozzles  were  within  a  few  nozzle  widths  of  each  other.   For  the  same 
reason  the  nozzles  were  also  required  to  be  within  a  few  nozzle  widths 
of  the  ground.   The  nozzles  were  assumed  to  be  perpendicular  to  the 
ground.   Since  tilting  of  nozzles  is  one  possible  solution  to  recircu- 
lation problems,  a  more  general  analysis  should  allow  the  nozzles  to  be 
at  some  angle  to  the  perpendicular. 

Considering  the  flow  field  shown  in  Figure  2  and  proceeding  exactly 
as  did  Goldstein  and  Siegel  the  hodograph  and  complex  potential  planes 
are  as  shown  in  Figure  3-1  end  3-2.   The  function  which  maps  the 
rectangular  region  T   in  the  T-plane  (figure  3-3)  into  the  hodograph 
plane  is 

G(T)  =  Vo  ei(TT/2'a)  n(T;§H,k)  VQ(T;§G,k)  [0(T;§A>k)  ]3/n       (1) 

where  the  Q  function  is  w  defined  in  Appendix  A.   Taking  the  argument 
of  this  function  and  requiring  that  the  direction  of  the    flow  along  the 
plate  must  equal  the  argument,  it  can  be  shown  that 
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(1)   Hodograph   plane,    Q   =  u-iv 
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(2)   Complex   potential   plane,   W  =  $+ii|i 
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(3)   T-plane 
FIGURE  3 
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The  remainder  of  the  analysis  proceeds  exactly  as  before  and  the 
final  solution  for  displacement  in  the  flow  field  becomes 


z   /'Vr\   ia  k' 
A  =  K5FJ   e   " 


sn  (TiF,k')  +  6R/6L-, 


1  +  6D/6T 
o  R   L 

T  (3) 

(k  sn  §u  snT-1)  (1-k  snT)  dT 
X 


dnTll-dnCn  ,k')  snT]  Q(T;?  k)  VfXT;Sr,k)  [Q(T;?  k)]3/lT 
§„+iK'  F  H  G        A 

G 

There  are  in  this  case  six  parameters  which  determine  the  flow  field: 
plate  height,  plate  width,  plate  angle,  nozzle  height,  nozzle  spacing  and 
nozzle  angle. 

The  plate  and  nozzle  angles  appear  in  the  solution  but  t  2  other 
parameters  must  be  calculated.   Following  the  exanple  of  Goldstein  and 
Siegel  the  solution  is  determined  by  integrating  along  the  sides  of  the 
rectangular  region  T   in  the  counterclockwise  direction.   Thus  the  distance 
between  two  points  in  the  physical  plane  is  the  difference  between  the 
line  integral  from  the  point  G  in  the  T-plane  (origin  in  the  physical 
plane)  to  each  of  the  points  in  question.   If  all  distances  are  non- 
dimensionalized  by  the  nozzle  width  the  solution  can  be  obtained  using 
the  known  plate  and  nozzle  angles  and  the  four  following  input  parameters: 
k,  6T/6R,  rjp/K1  .  and  5»/K.   These  parameters  all  have  values  ranging 
between  zero  and  one,  but  no  representative  values  were  given  in  Ref.  3. 
Upon  direct  communication  with  the  Lewis  Research  Center  it  was  learned 
that  no  values  of  these  parameters  and  their  corresponding  physical  con- 
figurations were  available. 
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IV.   COMPUTER  ANALYSIS 

Using  the  working  formulas  in  the  appendixes  of  Refs .  3  and  4,  two 
computer  programs  were  written.   The  first  program  calculated  the  plate 
height  and  width  and  nozzle  height  and  spacing  for  a  given  set  of  input 
parameters.   The  second  program  was  merely  a  variation  and  extension  of 
the  first.   For  a  given  set  of  input  parameters  and  experimental  data  the 
second  program  plotted  the  theoretical  flow  field,  the  pressure  coeffici- 
ents on  the  ground  and  plate  and  the  velocity  profile  in  the  nozzles  along 
with  the  data  points.   A  copy  of  the  second  program  follows  Appendix  A. 

The  plate  and  nozzle  angles  were  read  from  data  cards.   In  the  second 

program  the  input  parameters  were  also  read  from  cards,  but  in  the  first 

program  the  parameters  were  initialized  and  iterated  in  steps  of  one 

tenth.   The  complete  elliptic  integral  of  the  first  kind,  the  complimentary 

modulus  and  the  complimentary  integral  were  all  calculated  from  k.   The 

various  quantities  and  functions  were  then  evaluated  as  indicated  in  the 

program. 

The  modified  equation  for  §   is  equation  (2)  of  this  thesis.   The 

G 

modified   equation    for  0   is 

®(i)(^)  =f  -  a+  Q^^T!)  +|Q^i)(T1)  +^Qf)(r1)     1-3,4         (4) 

These  modifications  provide  for  solutions  with  the  nozzles  tilted.   The 
Q    's  are  angles  defined  in  Ref.  4.   It  should  be  pointed  out  that  Q 

must  be  taken  as  a  negative  angle  in  the  third  or  fourth  quadrant,  while 

(4) 
Q    must  be  taken  as  a  positive  angle  in  the  first  or  second  quadrant. 

The  height  and  width  of  the  plate  were  calculated  directly  by  inte- 
gration.  The  no:.zle  spacing  was  found  by  first  locating  the  poir,  B  at 
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the  end  of  the  plate  and  then  adding  to  that  the  horizontal  separation  of 
the  points  B  and  C.   The  nozzle  height  was  found  similarly  by  adding  the 
vertical  separation  of  the  points  D  and  E  to  the  width  of  the  flow  to  the 
right . 

At  this  point  the  values  calculated  were  printed  out.   In  the  first 
program  the  parameters  were  iterated  here  and  a  new  configuration  was 
calculated.   In  the  second  program  plotting  was  begun  at  this  point. 

The  data  points  were  plotted  first  and  then  the  theoretical  free 
streamlines  were  plotted.   A  similar  procedure  was  followed  for  the 
pressure  coefficients  and  the  nozzle  exit  velocity.   The  velocity  was 
made  negative  to  fix  the  direction  downward,  so  that  the  plotted  profile 
was  as  if  it  were  being  viewed  from  in  front  of  the  aircraft. 

All  integrations  were  performed  by  the  Gauss  nine-point  quadrature 
formula.   When  infinite  series  were  specified  only  the  first  five  terms 
were  taken.   In  the  working  formulas  presented  in  Ref.  3  quite  often 
square  roots  and  odd  powers  were  taken  of  the  negatives  of  numbers.   On 
the  intervals  indicated  this  would  be  quite  proper  since  the  numbers 
would  be  negative  values  before  having  the  sign  changed.   Since  the 
programs  were  to  be  used  to  find  possible  combinations  of  the  input 
parameters,  quite  often,  if  the  parameters  were  wrong,  the  result  would 
be  an  attempt  to  take  the  logarithm  or  square  root  of  a  negative  number. 
To  prevent  the  computer  from  ending  the  calculation  with  error  messages, 
absolute  values  were  taken  of  all  functions  under  radicals.   This  was 
equivalent  to  the  operation  intended  in  the  original  notation,  but 
avoided  the  problems  in  calculation. 

In  the  first  program  statements  were  included  in  all  the  functions 
to  notify  the  operator  if  a  denominator  was  zero.   The  second  program 
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has  retained  the  check  on  the  denominator  and  sets  the  function  to  a 
nominal  value  of  one  thousand  in  the  event  the  denominators  are  very 
small.   This  is  intended  only  to  avoid  underflow  errors  at  the  limits  of 
the  integration;  no  notification  is  given  to  the  operator  in  this  case. 
The  second  program  is  therefore  intended  to  be  used  only  with  combinations 
of  the  parameters  which  are  known  to  satisfy  all  the  equations. 
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V.   DESIGN  CONSIDERATIONS 

A  model  of  a  V/STOL  aircraft  using  two-dimensional  slot  jets  in  an 
inviscid,  incompressible  fluid  was  required  to  verify  the  theoretical 
results.   The  test  apparatus  was  designed  to  approximate  these  conditions 
as  closely  as  possible. 

Air  at  low  pressure  was  used  to  represent  the  incompressible,  inviscid 
fluid.   The  two-dimensional  slot  jets  were  produced  by  two  rectangular 
nozzles  with  a  cross -sectional  width  to  length  ratio  of  four  to  one.      By 
confining  all  measurements  to  a  plane  through  the  centers  of  the  two 
nozzles,  end  effects  could  be  ignored. 

To  avoid  entrainment  problems  the  theory  required  that  the  nozzle 
half-spacing  and  height  be  only  a  few  jet  widths.   This  condition 
suggested  the  use  of  a  large  nozzle  width  so  that  there  would  be  enough 
room  to  probe  the  flow  field  beneath  and  between  the  nozzles.   Also  to 
permit  a  significant  number  of  data  points  for  plates  with  half  widths 
as  small  as  one  fourth  of  a  jet  width,  a  large  nozzle  width  was  required. 
The  test  apparatus  was  designed  with  nozzles  two  inches  wide. 

Even  with  low  nozzle  exit  velocities  a  larg;  flow  rate  was  necessary 
to  supply  the  system.   Thus,  the  nozzles  were  connected  to  the  laboratory 
turbine  compressor.   Each  nozzle  was  attached  to  the  air  supply  through  a 
flow  straightening,  diffuser  chamber.   Rotameters  preceded  these  chambers 
in  the  air  lines.   Thus,  air  flow  conditions  could  be  monitored  easily 
throughout  the  test  period.   As  air  entered  the  diffuser  chamber  it  was 
deflected  by  three  layers  of  wire  screen.   The  air  then  passed  through  a 
layer  of  honeycomb  material  into  the  nozzle,  thus  producing  a  uniform 
velocity  distribution  throughout  the  nozzle. 
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Laboratory  model   of  V/STOL  aircraft   using 
symmetrical   two-dimensional   slot   jets 

FIGURE  4 
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Flow  straightening  chamber 
FIGURE   5 
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The  flow  straightening  sections  also  served  as  a  mounting  for  the 
nozzles.   If  desired,  the  nozzle  shape  and  size  could  be  changed  quite 
easily . 

A  mounting  arrangement  was  designed  to  not  only  allow  the  nozzles  to 
be  moved  vertically  and  horizontally  but  also  to  be  rotated  in  the  vertical 
plane,  thereby  representing  a  tilted  jet  configuration.   This  freedom  of 
movement  was  accomplished  by  connecting  the  flow  straightening  section  to 
a  vertical  rod  through  a  swivel  which  could  be  locked  at  any  angle.   The 
vertical  rod  passed  freely  through  a  metal  block,  and  it  too  could  be 
locked  in  any  position.   The  block  moved  horizontally  along  parallel  bars 
and  could  be  held  firmly  in  position  with  a  set  screw.   The  parallel  bars 
were  mounted  in  a  rigid  framework  which  was  supported  by  four  legs. 

Directly  underneath  the  framework  and  the  nozzles  was  a  large  plexi- 
glass plate  which  represented  the  ground.   Smaller  plates  representing 
fuselages  of  various  sizes  and  shapes  were  made  to  be  interchangeable  on 
a  bracket  which  was  suspended  from  the  framework  between  the  nozzles. 
The  vertical  position  of  the  plate  bracket  was  varied  with  a  threaded 
traverse.   The  ground  plate  and  the  fuselage  plates  contained  a  row  of 
pressure  taps  across  the  center  section.   For  future  studies  with  smaller 
nozzles  and  higher  velocities,  significant  pressure  measurements  would  be 
made  possible,  but  for  these  low  velocity  tests  the  primary  measurements 
were  made  with  a  hot-wire  anemometer. 

To  support  the  hot-wire  probe  a  traversing  mechanism  was  mounted  on 
the  apparatus.   The  mechanism  allowed  the  probe  to  be  held  at  a  constant 
depth  into  the  flow  while  traveling  both  horizontally  and  vertically 
throughout  the  field. 
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Nozzle  mounting  arrangement 
FIGURE  7 


Traversing  mechanism 
FIGURE  8 
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VI.   PROCEDURE 

A  list  of  plate  heights  and  widths  and  nozzle  heights  and  spacings 
for  which  the  corresponding  input  parameters  were  known  was  compiled. 
From  this  list  a  few  representative  configurations  were  chosen  for 
experimental  verification. 

The  plate  size,  height  and  angle  and  nozzle  spacing,  height  and  angle 
were  all  set  and  recorded.   The  compressor  was  then  started  and  an  equal 
air  flow  was  routed  to  each  nozzle  through  the  flow  meters.   A  straight 
hot-wire  probe  was  placed  under  one  edge  of  one  of  the  nozzles.   The 
voltage  and  location  was  recorded.   Then  by  following  the  motion  of  a 
fine  thread  placed  in  the  flow  field,  the  probe  was  moved  horizontally 
in  convenient  intervals  and  then  moved  vertically  until  the  same  voltage 
was  again  observed.   Since  only  velocity  magnitudes  were  of  interest,  the 
probe  was  always  rotated  to  obtain  a  maximum  reading.   By  following  this 
procedure  from  both  sides  of  a  nozzle  and  from  the  edge  of  the  upper  plate 
the  coordinates  of  the  lines  of  constant  voltage  (free  streamlines)  could 
be  recorded. 

The  probe  was  then  moved  as  close  to  the  upper  plate  as  possible  and 
oriented  to  read  the  velocity  in  the  direction  of  the  plate.   By  moving 
the  probe- in  convenient  intervals  along  the  plate  the  velocity  distri- 
bution along  the  plate  couM  be  recorded.   Similarily  the  probe  was  moved 
along  the  lower  plate  and  the  voltage  and  horizontal  position  were 
recorded . 

This  procedure  was  followed  with  each  configuration.   Tests  were  run 
with  Reynolds  numbers  varying  from  1600  to  7000.   Simultaneously  pressure 
measurements  were  made  using  a  micromanometer  capable  of  indicating 
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pressures  as  low  as  .001  inches  of  water.   The  pressure  and  location  were 
recorded  for  comparison  with  the  hot-wire  data. 

The  data  collected  in  this  manner  was  sufficient  to  plot  the  free 
streamlines  of  the  flow  field,  the  nozzle  exit  velocity  profile  and  the 
pressure  coefficient  along  the  plate  and  ground.   The  input  parameters  for 
each  configuration  were  read  into  the  program  along  with  the  experimental 
data.   The  experimental  results  for  each  configuration  were  then  plotted 
along  with  the  corresponding  theoretical  results. 
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VII.   RESULTS 

The  results  of  the  experiments  are  presented  in  the  following  pages. 
Table  I  is  a  list  of  some  of  the  input  parameters  for  which  the  corres- 
ponding configurations  were  found. 

Figures  8-1  through  12-1  are  typical  free  streamline  plots  for  various 
configurations.   Hot-wire  data  points  are  represented  by  triangles.   The 
coordinates  are  nondimensionalized  with  the  nozzle  width. 

The  pressure  co  Cficient  along  the  ground  is  shown  in  Figures  8-2 
to  12-2.   Only  the  abscissa  is  nondimensionalized.   The  pluses  represent 
data  points  calculated  from  manometer  readings. 

The  pressure  coefficient  on  the  plate  is  shown  in  Figures  8-3  through 
12-3.   Again,  only  the  absciss,  is  nondimensionalized. 

Figures  8-4  through  12-4  are  the  plots  of  the  nozzle  exit  velocity 
profile.   The  velocity  is  nondimensionalized  with  respect  to  the  maximum 
velocity  recorded  for  that  test.   The  abscissa  is  one  nozzle  width  in 
length.   The  velocity  is  plotted  negatively  so  that  the  plot  represents 
the  actual  direction  of  the  flow. 

Figures  13  and  14  clearly  indicate  the  nature  of  the  flow  field 
beneath  the  nozzles. 
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VIII.   DISCUSSION  OF  RESULTS 

Many  combinations  of  the  input  parameters  produced  impossible  con- 
figurations (negative  plate  heights,  etc.)-   In  fact  there  were  far  more 
extraneous  solutions  than  useful  ones.   There  were  no  combinations  found 
which  produced  a  configuration  wi  h  the  nozzle  exit  lower  than  the  plate. 
Very  few  reasonable  configurations  resulted  from  the  maximum  or  minimum 
values  for  the  parameters . 

Most  of  the  configurations  of  interest  were  found  to  have  input 
parameters  in  the  following  ranges : 

.3  <  k  <  .6 
.1<6L/6R<  .5 
•6  <  Tip/K'  <  .9 
•1  <  ?A/K<  .9 

There  were  possible  configurations  outside  of  these  ranges,  but  most  of 
them  had  extremely  small  plate  widths.   Variation  of  any  one  input 
parameter  was  not  found  to  change  any  single  physical  parameter. 

Decreasing  the  plate  angle  had  a  tendency  to  decrease  the  plate  width 
and  increase  the  plate  height.   Tilting  the  nozzles  caused  an  increase  in 
plate  and  nozzle  height  but  decreased  plate  width  and  nozzle  spacing  for 
a  given  set  of  abstract  input  parameters.   Neither  of  these  variations  had 
any  outstanding  effect  on  the  theoretical  pressure  coefficients. 

The  experimental  data  for  the  streamlines  agree  quite  well  with  the 
theoretical  plot.   The  data  points  between  the  plate  and  nozzle  often 
seem  to  follow  a  different  streamline  than  the  theory.   This  is  a  result 
of  the  fact  that  the  test  plates  were  of  finite  thickness,  while  the 
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theoretical  plates  were  infinitely  thin.   Figure  9-1  shows  entrainment 
effects  for  very  high  nozzle  positions.   In  some  cases  the  theoretical 
plot  can  be  very  impractical.   Notice  in  Figure  11-1  that  the  streamline 
from  the  edge  of  the  plate  turns  and  passes  through  the  plate. 

In  most  cases  the  theoretical  pressure  coefficients  along  the  ground 
and  plate  were  quite  close  to  the  actual  situation.   The  micromanometer 
data  points  are  somewhat  questionable.   The  pressures  were  extremely  low 
even  for  the  instrument  which  was  used.   The  theoretical  coefficient  is 
assumed  to  approach  zero  after  the  last  point  plotted,  as  the  measurements 
indicate . 

The  nozzle  exit  velocity  profiles  for  very  low  velocities  were 
approximately  the  same  shape  as  those  from  the  theory,  but  as  the  exit 
velocities  become  higher  the  profile  becomes  much  flatter.   In  all  cases 
it  should  be  remembered  that  the  velocity  data  points  are  normalized  with 
the  maximum  velocity  for  that  particular  test.   In  most  cases  this  maximum 
velocity  was  found  to  be  near  the  ground  two  or  three  nozzle  widths  from 
the  high  pressure  region  beneath  the  fuselage  and  nozzle  exit.   Hence,  the 
data  points  tend  to  be  displaced  from  the  theoretical  profile  which  assumes 
that  the  maximum  velocity  is  reached  at  the  edges  of  the  nozzle  exit.   In 
all  cases  the  profile  was  somewhat  affected  by  the  presence  of  the  wall  of 
the  nozzle. 

The  motion  of  the  tufts  of  yarn  placed  in  the  nozzle  exit  (Figure  13) 
show  the  manner  in  which  the  flow  divides  to  the  left  and  right  of  the 
nozzle.   The  strips  of  paper  attached  to  the  nozzle  exit  and  ground  plate 
in  Figure  14  show  the  stream  of  flow  away  from  the  nozzle  an   the  upward 
air  currents  under  the  fuselage. 
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IX .   CONCLUSIONS 

The  experiments  have  not  only  verified  the  theory,  but  have  also 
pointed  out  its  limitations.   The  analytical  results  were  quite  accurate 
for  nozzles  very  close  to  the  plate  and  the  ground.   As  the  nozzles  were 
moved  higher  the  effects  of  entrainment  became  more  noticeable.   The 
pressure  distribution  along  the  ground  was  not  as  uniform  as  that  pre- 
dicted by  the  analysis. 

The  greatest  limitation  on  the  theory,  however,  is  not  the  restricted 
nozzle  positions,  but  the  fact  that  the  analytical  solution  depends  upon 
four  parameters  in  an  abstract  plane.   Since  there  is  no  direct  connection 
between  any  one  of  these  pai'ameters  and  the  physical  parameters  of  plate 
height  and  width  and  nozzle  height  and  spacing,  the  theoretical  results 
are  only  possible  after  an  iteration  process  to  find  the  proper  input 
parameters.   A  vast  amount  of  computer  time  is  required  to  solve  for 
relatively  few  configurations. 
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X.   RECOMMENDATIONS 

The  theoretical  solution  could  possibly  be  extended  to  include  other 
fuselage  shapes.   However,  perhaps  a  more  practical  approach  would  be  to 
develop  analytical  solutions  in  terms  of  physical  parameters.   A  solution 
of  this  type  would  be  a  useful  tool  to  the  aircraft  designer. 

Further  tests  should  be  carried  out  with  higher  nozzle  exit  velocities 
to  determine  the  effects  of  extremely  turbulent  flow  underneath  the  air- 
craft.  Also  of  interest  would  be  the  effect  of  crosswinds  on  the  flow  and 
the  pressure  distribution. 

Reference  1  mentions  the  use  of  a  special  takeoff  grating  to  virtually 

eliminate  entrainment  effects  on  aircraft  using  a  single  central  jet 

exhaust.   No  information  was  available  concerning  the  effects  of  such  a 

grating  on  multi-jet  configurations.   Tests  of  a  takeoff  grating  mounted 

on  the  existing  apparatus  would  be  of  interest  to  determine  the  resulting 

changes  in  ground  effect. 

There  were  no  reasons  given  for  restricting  6,/&D  to  values  less  than 

L      K 

one.      Obviously  there  must   be   solutions    for   nozzle  exits    lower   than   the 

plate   height.      If  the   equations   can  be   satisfied   by  values    of  6,/&D 

L  is. 

greater  than  one,  perhaps  this  will  lead  to  the  lower  nozzle  positions. 
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APPENDIX  A 

To  obtain  analytical  solutions  for  the  flow  field  beneath  a  V/STOL 
aircraft  Siegel  and  Goldstein  [Ref.  3]  considered  a  simplified  two- 
dimensional  model  of  the  flow.   The  following  is  a  condensation  of  their 

analysis  . 

The  aircraft  under  consideration  has  a  fan  pod  or  jet  exhaust  in  each 
wing.   As  the  two  jets  strike  the  ground  they  divide  into  flows  to  the 
left  and  right.   The  colliding  flows  beneath  the  aircraft  rise  and 
strike  the  fuselage.   The  bottom  of  the  fuselage  is  represented  by  a 
flat  plate  which  can  be  at  an  angle  with  the  ground.   Because  of  symmetry 
only  half  of  the  flow  field  is  considered. 

The  flow  is  assumed  to  be  incompressible,  inviscid,  irrotational , 
isothermal  and  steady.   The  velocity  at  the  nozzle  exit  is  taken  to  be 
perpendicular  to  the  exit  plane.   The  surrounding  fluid  is  at  rest  and 
entrapment  is  ignored.   If  the  nozzle  exits  are  within  a  few  jet  widths 
of  each  other  and  of  the  ground,  the  inviscid  solution  with  no  entrap- 
ment is  a  reasonable  approximation  of  the  actual  flow  conditions.   With 
these  restrictions  the  flow  has  been  reduced  to  a  free  streamline  prob- 
lem with  constant  pressure  surroundings. 

Figure  Al-1  shows  the  physical  boundaries  of  the  flow  field.   Defining 
the  complex  conjugate  velocity  as 

C  =  u  -  iv  <A1> 

and  the  complex  potential  as 

W  =  $  +  if  (A2) 

the  flow  regions  in  the  hodograph  and  complex  potential  planes  are  shown 
in  Figures  Al-2  and  Al-3. 
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The  solution  in  the  physical  plane  is  obtained  by  performing  the 
integration 

z  =  /  —  ciW  +  constant  (A3) 

In  order  to  perform  this  integration  Q   and  W  must  be  expressed  in 
terms  of  the  same  variables.   This  is  accomplished  by  mapping  the 
rectangular  region  T   of  the  T-plane  onto  the  Q   and  W  planes. 

The  function  which  maps  the  region  T   onto  the  hodograph  plane  is 

C(T)  =  Voei(!T/2)Q(T;§H,k)  VQ(T;§G,k)  [fi(T;§A>k)  ]p/n       (A4) 
The  Omega  functions  are  defined  in  Ref.  4  to  be 

Q(T;§  ,k)  =  — j- — -FT  (A5) 

6Sm  <t+V2k>I§-) 

where  9,  is  the  fourth  theta  function  which  is  defined  in  Ref.  8  as 
4 

follows : 

9,  (z|t)  =  1  +  2   E  (-l)n  eim  cos  2nz  (A6) 

n=l 

The  relationship  between  §.  ,  £,      and  §  must  satisfy  the  requirement 

AG       H 

that  on  the  plate  the  argument  of  £(T)  must  equal  the  direction  along  the 
plate.   By  using  the  formulas  for  the  arguments  of  the  Omega  functions 
presented  in  Ref.  4  it  can  be  shown  that 

By  applying  a  Schw..:rz-Christof fe  1  transformation  the  upper  half  of 
an  intermediate  t-plane  (not  illustrated)  is  mapped  onto  the  complex 
potential  plane.   A  second  transformation  maps  the  T-plane  onto  the 
upper  half  of  the  t-plane.   By  eliminating  t  between  these  two  trans- 
formations, the  function  which  maps  the  T-plane  onto  the  W-plane  is 
found  to  be 
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dW    ..     (k  Sn  ^H  Sn  T  -  *>  dn  T 

dT  "'  lA  (1  +  k  sn  T)  [l-dn(TU,k!)  sn  T]  (A8) 

r 

where  A  is  a  constant. 

Using  partial  fractions,  integrating  equation  (A8)  and  applying  the 
condition  that  the  flow  through  the  nozzle  exit  is  equal  to  the  differ- 
ence in  the  stream  function  across  the  exit,  the  value  of  A  can  be  found 
in  terms  of  the  average  nozzle  exit  velocity. 

Noting  that  the  flow  to  the  left  beneath  the  nozzle  is  equal  to  the 
change  in  the  imaginary  part  of  the  complex  potential  at  F  and  requiring 
the  velocity  at  that  point  to  be  the  free  streamline  vc  locity  we  have 

k  +  dn  (T]F,k')  -  k  (1  +  sn  ?H) 
6LVo  =  V  k  +  dn  (T]F,k')  -  k  (1  +  sn  §H)  [1  -  sn  (^.k1)]   (A9) 

Since  continuity  requires  that 

6nV  +  6TV  =  AVT  (A10) 

R  o     L  o     I 

we  can  eliminate  6,  from  equation  (A9)  and  have 

Ju 


6R        V_  k    (1  +  sn  F    )    sn    (t^k') 


R  I  3H 


A         Vo  k  +  dn    OlF,k')    -   k   (1  +  sn  §R)    [1   -   sn    (r,F>k')] 
By  dividing  equation    (All)   by  equation    (A9)   and  rearranging 


(All) 


6R/6T    dn    Oi^k')    -   k  sn    (TU,k») 


R'     L 
iR/tL 


Sn  ?H  =        6D/6T    k  +  k   sn    (r,      k')  <A12> 


Since   Q   and  W  are   both  known   in   the   T-plane   the   solution   of  equation 

(A3)    becomes 

T 

■-  f     tkf?dT  '  (A13) 

5G+iK' 
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or  after  substitution 


V. 


:n  (nF,k«)  +  6R/6L- 


i>  hi 

A  "  VV  J   tt  L    1  +  6„/6 


R'  L 


(k  sn  5R  sn  T  -  1)  (1  -  k  sn  T)  dT 


(A14) 


?G+iK' 


dn  T  [l-dn(nF,k')  snT]  Q(T ^R,k^^f7^^    [Q(T;?A,k)] 


3/tt 


Using  this  relation  the  ratio  of  the  average  velocity  in  the  nozzle 
exit  to  the  free  stream  velocity  is 


V 
o 


1  +  V6L 


Lsn  (TiF>k')  +  6R/6L- 


(A15) 


where 


K 


M  =  < — 


(1-k  sn  ?H  snF)  (l-k  sn§)  d^ 


{n  j  dn5[l-dn(nF,k')sn5]n(§;§H,k)Vn(§;§Gjk)[Q(§;?AJk)]P/n 
-K 


(A16) 


The  pressure  coefficient  at  a  point  (x] ,  y  )  is 


C   = 


p<v  y{>  -  p0 


P    1/2  p  V 


(A  17) 


Using  Bernoulli's  equation  the  pressure  coefficients  on  the  ground 
and  on  the  fuselage  are 

P,8 


c       -  i  -  [H-aJQl  -k  <  §  <  5 

p,g  L       V         J  -  * 


and 


(A18) 


C      .  =  1 
P.f 


[iL-Jg.K^Y    <MO]       5a  <  5  < 


<  K 


The  flow  field  is  governed  by  five  physical  parameters :   plate  height, 
plate  width,  plate  angle,  nozzle  height  and  nozzle  spacing.   Only  the 
plate  angle  (3  appears  in  the  solution.   Hence  the  other  parameters  must 
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be  found  in  terms  of  four  parameters  in  the  T-plane.   Siege  1  and  Goldstein 
used  k,  6j/6R,  tu/K'  and  £A/K  as  input  parameters  for  a  computer  program 
to  plot  numc  Leal  results.   All  four  parameters  have  values  between  zero 
and  one  but  no  representative  values  and  their  corresponding  physical 
configurations  were  available. 
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